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Abstract 


A plasma simulation code is presented that models the time-dependent plasma 
properties in the vicinity of a Spherical, charged spacecraft. After showing 
agreement with analytic, steady-state theories artd ATS-6 satellite data, the follow- 
ing three problems are treated: (1) transient pulSes from photoemlSslOn at various 
emission temperatures and ambient plasma conditions. (2) spacecharge limited 
emission, ^ and (3) simulated plasma oscillations in the long- wavelength limit 

(k kjj « 1). 



1. INTRODUCTION 
1.1 Mbtfel Objecllvbs 

The objective of this computer model is to realistically simulate plasma- 
spaceCraft interactions. It presently treats time-dependent plasma phenomena ifi 
the limit of spherical symmetry. Although future plans anticipate the incorporation 
of a realistic three-dimensional spacecraft geOmetry, Understanding Of the spheri- 
cally Symmetric limit uniquely identifies plaSma effects. In this manner, it is 
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ht>i>ed that the wodel will pi'ove to be a uaeful tool Itl differentiating between plasifla 
and geometric phenomena observed in the SCATHA Satellite data. 


1.2 OoHcrlpU^lt of iho Mdddl 

We use what is commonly called a "p&«lcle pusher" model, that Is, one 
approximates the actual plasma by a number of "computer" particles whose charge, 
position, angular momentum, ana velocity are tracked in time. By properly 
weighting these computer particles, reasonable statistics can be Obtained near the — 

spacecraft. ^ 

Appropriate particle aistrlbutions are generated either by a Monte Carlo 
technique or by a systematic loading of the velocity and spatial intervale ("quiet 
start"). ^ 

For purposes of comparison, we have initially restricted ourselves to 
Maxwellian distributions. However, both methods (Monte Carlo and Quiet Start) 
can easily be extended to any distribution that can be numerically integrated. Once 
the neutral plasma IS created, the Computer tracks the particles and recalculates 
the potential at each time Step. Particles that hit or are emitted from the space- 
craft are taken Into account as well as those that enter and exit the sheath boundary. 
Thus, tiie eomputeriaed sheath structure dynamically eVolVeS in analogy with the 
physical Situation. Ohe may then store the results and restart the program With a 

new set of environmental parameters. 

In Section 2 of this paper, we will show the agreement of the present work 
with that of other approaches and with the ATS- 6 data. Section 3 deals with the 
photosheath. The minimum rise time of the Spacecraft potential Is determined for 
various emission temperatures and ambient plasma parameters. Space charge 
effects are alSo discussed as Well aS those due to Secondary emission and back- 
scattering. Finally, in Section 4, we Illustrate the possible presence of plasma 
oscillations by performing "computer experiments" with the Code. 


2. comparison with analVtic approaches 

2.1 CotDparlmii with Latigmtilr Theory 

As a first step We compared our steady-state results with those expected 
from Langmuir theory.^* ® figure 1 shows the comparison, in these results, the 
ambient plaSmS temperature and density were held fixed and the fixed potential on 
the probe (^acecraft) Increased. The parameter Ip Is the ambient current to the 
probe surface at aerO Voltage, if one makes the Same assumption fbr Oie mOdel aS 
for the Langmuir theory (that Is, rtO presheath acceleration), good agreement IS 
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Xq iS-SSm 


MOtT'SMlTH > LANOMUIA THEORY 
• •SIMULATION RROORAM WITH RRESHEATH ACCELERATION 
X - SIMULATION RROORAM WITHOUT RRESHEATH ACCELERATION 



Figure 1. Cdnfiparlsbii of the Sithuietlon Code with the Mott>Smith- 
Langitiuir Theory. This figure depicts a current-vOltage cUrve iil 
the thick'Sheath approximation, l^e probe radius is 1 Ri and the 
outer sheath boundary has been set to Fq = 2 and = 3 in units of 
k] 3 . The dots represent the ceae where there is presheath acceler- 
ation. That IS. when the potential at Rg is nonzero 


Obtained. (Presheath acceleration takes into account the dependence of the sheath 
size on the probe potential. In the code, this is represented by a Roltzmann 
factor at the sheath boundary. 
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2.2 

we eleo eompered our reeultd 1.. ate.dy-tBto lle,lt «, „unerl<-al »olu«one 

Jvlovadd Peleeonofld.tiene.'’'’ The d.e»d, UmiU. re.edod d. 

Figure 2 dlu.»e the code reoulta alter 0. 12 m.oe. In thl« cample, th. "» 

a bU9 rollngc equal to mlnua ten times the a.erage electron energy cpreae- d m 

oteetron ,olto. From Figure 2 it l. eeen ZiZZ°'nr‘ 

prollle. but dlHor eome»hat lor the electron denelty nt large r Tl.t. dill 
L probably due to the lime-dependent code not having reached the true eteady 

’““FCrVa aho.e . elmllar eomparleon lor relatively intenee monoenergetle 
etee^emleelon Irom the probe surlace. « The probe la poaltlvel, btaaed at 
« 0 v^ The emlaaton energy la *1. 0 electron volta ao that at. 
trona return to the probe auMace. The re.ultant denelty 

adreetnent. Scatter In the time -simulation results can be Improved by using 
Lre computer partlclea. For cample, tn the preaent case «e >■"'* =‘PP~»""“ ' 
rooo computer tone and eleetrona. Thin number can be aignllteantly u..-rea.ed at 

the expense ot longer computation time. 

2.3 Coirtparlson with ,\TS-6 Daia 

TWO charging events in the eclipse region wfere analyzed from ATS-6 data 
aupplled by DeForeat. » The proton and electron ^ 

as well as the simultaneous vi'hicle potential are nhown i , ,, • 

evlto rireaent a time .hen Iniectton ol hot plaama tool, place. Alao aho.n m 
Table 1 aM the vehicle potenllala as predicted Irom a almple Boltamann approx, 
on «d a toleh-aheath approxl. a«on to the orblt-llmlted Langmuir theory. 
The Boltamann appro- Im.tton aaaumea that the electron denelty .a ctoaely 

repr.aentedb,N,-N^e-'*P' ' Protona are conaldered to be unatlected. 

Therefore, current balance is given by 


kT - I eeJ /kT 






at 


eyktg ■= tnUtgfnp/tpm^l 


1/2 


( 1 ) 


where mp = protort mass, m^ 
defined below. 


electron mass, k = Boltamann constant and c is 
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— PARKER'S ■PARKLAP*- R,»3 
“.AAFSL'S "SHEATH* R. >3 
(t«l,TKlO"*»#e) 



Figure 2. Potential and Density Profiles in the Sheath. Com- 
parison of the time-dependent cdde with a Steady-state numeri- 
cal solutidn to the Vlasov-PoiSSon equations (L. W. Parker, 
private communieation). The density, n, is normaliaed to its 
ambient Value. The potential, V, is showh norm&Uaed tb the 
electron temperature. The parameter, t , is the time at 
which the time-dependent results were taken. A negative bias 
of -10 kT- volts is on the Spacecraft SuMace, Rn is exptessed 
in units of probe radii 
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Table 1. Comperleon ot BollMoannjmd. Langmuir Thenrlea .lib ATS-8 
DAtS. € “ 0. 16 


Dfev of Event 
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(1976) 




e^l/kT. 


MeAS. Boltzmann 


23. "i 9*1 
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'the L^ngmuil* tWdk-fiheath apprdXlm&ttoki is more realietie in that it also 
takes Into account the effect of the Spacecraft potential on the net proton curreht. 
It Is given by 


-U^pl/kT 

\7v 



(1 + eyUTp) 


( 2 ) 


a trahScendental equation in4he eatellUe potential. A . which can4>e reduced to a 
quadratic equation in The parameter, e, ’repi*esentfi the net fractional 

electron current to the spacecraft and, thus, takes into account secondary emis- 
sion. backscattering and photoemissioh. In Table I we have set e = 0. 16 (that is. 
fi4 percent net backSCatter. etc. ) for both approximations. Examination of Eq. (2) 
shows that in the limit of large proton temperatures the thick-sheath approxima- 
tion reduces to the Boltzmann case. This is also seen from Table 1. In both 
cases, the Langmuir thick-sheath approximation gives better agreement with 
measured results than the Boltzmann limit. 

Insight into the expected sensitivity of spacecraft Uoltage to changes in the 
ambient current, can- be seen from Figure 4 which contains a plot of Eq. (2). The 
top-half of this figure is a linear Scale plot of normalized voltage to the net electron 
current incident on the spacecraft normalized to the ambient electron current 
(With the Boltzmann factor). The bottom half of the figure Shows the same curve 
On a log-log Scale plot. iTiese curves will be modified by any voltage dependence 
in c Such as Space charge limiting effects. 

In Conclusion, for large negative spacecraft potentials, the thick-sheath 
approximation predicts the average Spacecraft voltage, given the piroper material 
characteristics and the ambient electron and proton temperatures. The code is 
consistent with the thlCk-sheath limit, and predicts sheath density and voltage 
profiles. The simulation code, however, is also valid Where the thick-sheath 
approximation breaks down. That is, where space charge, and time -dependent 
effects become important. 

The importance of space charge is determined by Comparing the spacecraft 
surface charge to the charge residing in the sheath. For large surface potentials 
and v .nuous plasmas, spacecraft surface charge is dominant. In that case, 
Laplacian solutions with appropriate geometric boundary conditions should be 
adequate, At »oW vehicle potentials with surface emission, space charge effects 
become important. 
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THlCk - sheath APPROXlMAtiON Tp~T, 
linear PLC)f 



Ftgure 4 . Spacecraft Vpitage versus liet 

NprmSlized to the Ambient Current (Including the Boltzmann 

Factor) 


3. THE PHOfOSHEAtH 
3.1 InlroducUon 

incident Sunlight cauSes the etolsBlon of lew energy (~2 eV) electrons front 
the Sdtelllte surface, this emlsSldti. in a hot plasiha enyirdnirtCrit. soirtetlmes 
swlnas the sateliltfe botentlal thousands of volts as the satellite enters and emerges 





I t 


te<Jm th6 eclipse region. A second effect Of photoertitsfilon is Its background 
effect Oh ambient plasma density measurements. These measurements rnust be 
corrected for electrons originating on the satellite surface. A third effect of the 
photoSheath. Is its possible Interaction with onboard electron emitters either 
through space*charge limiting effects or by Stimulation of plasma oscillations. 

In the present code, we can simulate either mOnoenergetlc or Maxwellian 
photoemiSSion. The Maxwellian caSe is described in detail in Appendix A. The 
angular distribution of the emitted photoelectrons fbHuW a COSine law relative to 
the surface normal. In ttie following paragraphs, we consider photoemission 
to be switched "on" at t = 0. This is a worst case condition since satellites 
emerge from eclipse over a period of minutes. 

3.2 Monoenergetic i!mis.«inn 

Figure 5 Shows the simulation results for monoenergetic emission. The 
vertical axis represents the jatellite (taken as a 1 m radiUs Sphere) potential 
While the horizontal axis denotes time in miCrOSecOnds. The Satellite IS taken to 
be at zero volts at t = 0. The bump in each of the curves occurs when the initially 
emltted-electrons return to the satellite. The final surface is the result of two 
effects. First, emitted electrons do not return until the satellite reaches a voltage 
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equ&t t6 the ettitssloti-ehejigy. Secondly, thetb is & finite transit time for the 
returning electrons. This means that the final vdltage is Higher than the emission 
energy expressed in electron volts, the dotted lines represent Spacecraft voltage 
buildup in the limit Where all emitted electrons escape. Further cases regarding 
mohoenergetic photcemission can be found inJtata et al. 

3.3 tnittsie’iil-Bise*Tiiiie 

Figure 6 shows expected surface voltage rise -times at various ambient densi- 
ties. in these runs the photCelectrons are emitted with a Maxwellian distribution 
with a te*mt)6rature corresponding to 6 eV. The ambient defisity is then varied 
to detefctnlne the risfe-time^ senaitivity to the ratio of the photoelectrcn and arnbient 
currents. 


T«>O.S«V fephet >6«vnuiAXWeLLIAM) o • W ^ if.AMi 
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Figure 6. Spacecraft Voltage Transients fdr a Fixed Emission Temperature 
and-Various Flasma Densities. No secbndarieS inC'.uded 


In Figure t, we treat the alternative case. That is. the ambient plasma para* 
meters are held constant and the phdtdemiasion temperature is varied. The ddtted 
line represents 100 percent escape efficiency for the emitted electrons. The 
steady*8tate surface voltage, as expected is dependent on emission temperature. 

The transients shown in Figures 0 and "t Cbiild cause satellite malfunctions if 
they reached a critical logic circuit. Integrated Circuits usually Operate Over 
0*S vblts so that a ^ 2 volt transient through the grbund liiieS would give e faiSe 




* AMP/m* ij/cc T*^T( »d«V 



Figure 7. Transient Spacecraft Potentials for Strorlfr Photo- 
emission at V'^ arious Emission Temperatures. No secondaries 
included 


signal. The incident sunlight flux must change more rapidly than the characteristic 
time-constant of the surface materials. Otherwise, th^s^ rise times will not be 
significant Further research into this area needs to be performed. 


•i. 1 spare r.littrgi^ LiftiUed laiiissioh 

Chang and Bienkowski^^ showed that large current emission from a posi- 
tively biased prbbe is ihhibitbd by space -charge buildup in front of the surface. 

(Sed center Figure 8. ) Whipple^ Used U similar appi^oach to show that differential 
surface charging is necessary td explain the barrier potentials observed on AtS-f:. 
In this section, we use the simulatidn code to produce space-charge effects 
Application is then made to a floating spacecraft potential with and \^ilHdut active 
control experiments. 

The versatility df the present code is illustrated in the lefthand side of t'igure 8 
The surface potential is held fixed at +8 volts while the electrohdtnissidn current is 
increased. For a nominal 5 .10'^ A/m^ photoemisslon current density, no 
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potential bari*ier IS Obaefved. However, if the emission is artificially enhanced 
to 50 and SOOviO ’ A/m^ potential barriers are clearly evident and inhibit further 
emission, (in these runs, the sheath boundary was choSen at 3 32 rh) The effect 
of the potential barrier is sUcH that only 5t. 3 perOent. 23. 2 percent and 4 4 per- 
cent of the emitted flux reaches 3. 32 m for the three respective cases. These 
results Show that-the simulation code can treat space-charge effects and predict 
emission efficieneies. The code can distinguish between emitted and ambient 
particles and, therefore, can be used to correct ambient plasma measurements. 

The depth of the potential well is also dependent on the satellite voltage. On 
the righthknd side of Figure 8 we show three curves for the samt ambient and 
emission characteristics but at different surface potentials. For large negative 
Values Of the surface potential, the well is completely eliminated since emitted 
electrons are rapidly repelled. On the other hand, if the surface potential is too 
positive, the emitted electrons quickly return to the emitting surface and a signi- 
ficant Spacecharge cannot form. Therefore, potential wells are expected for only 
a Specific range of surface potentials. 

The situation iS more complejt with a floating spacecraft potential. In this 
case, both the well depth and the Surface Voltage are strongly dependent on the net 
electron current. Figure 9 shows the expected Voltage profile in the sheath at 
various emission currents dnd ambient densities. Runs were taken at three 
densities n = 1 Cm 5 cm ^ and 10 cm'^ with a photoemission current of 
8. 2p A/m^. The emitted current was then increased to 40 ^ A/m^ which is the 
limiting phOtOemiSsion intensity for aluminum. In addition to photoemission, 

84 percent of the incident electrons were assumed to cause iSotrOpic secondaries 
to be emitted with a Maxwellian energy distribution (kT = 2. .') eV). Tlie 84 percent 
figure was estimated from ATS-6 eclipse data (T&ble 1). Figure 9 shows that 
the presence Of a potential well is dependent On a "balance" between the ambient 
and emitted currents. Space charge also gives rise to the curious effect that a 
Spacecraft may have a net positive charge but be at a neg Uve potential relative to 
the ambient plasma. 

The Figure 9 also shows the maximum barrier to be about -3 volts. These 
results are in agreement With those of Whipple* who showed that the barrier 
potentials inferred from A IS -6 data are too large to be explained in terms Of a 
spherically symmetric photoelectron or secondary Sheath surrounding a uniformly 
charged spacecraft. Differential charging between spacecraft surfaces is, there- 
fore, mainly responsible for the potential barrier. The average satellite potential 
relative to the ambient plasma, however, is determined by the sheath. 

Active control experiments irt hot dense plasma could lead to more prOnoUrtced 
space-charge effects. As an extreme example, we took Tg = 9. 1 k eV, tj = 23.7k eV 
and n ^ 5v 10^ cm . In that case in order to maintain neutrality, electron emission 
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Figure 9. The Potential Prdfile in the Sheath at Various 
Emission Intensities and Ambient Densities. The sur- 
face potential Is floating 


of ~1. 6 A was reijuired. A potential barrier on the order of hundreds of volts 
was obtained. 

3..1 l‘'uiure 

The Monte -Carlo approach is presently being applied to include fealistir back 
scattering and secondary emission from both electrons and incident tons. These 
processes are energy and material dependent. In this way, insight into the impor- 
tance of rtiaterlal properties on spacecraft charging will be gained. 
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4. PU!^M\ OSClU.ATlOSS 


i. I — 4nUddi)ctioA- 

Because Of the dynarrtlc nature of the plasma, the Sheath potential may oscillate 
In time. Moving plasma particles "overshoot" theii> equilibritim posltlohs and 
uhdei*go simple harmonic motion. Collective oscillations Bite, therefore, closely 
connected with dynamic screening. The natural frequency of these collective 
modes is the plasma fjre<Juency, 



where n is the ambient plasma density. The sheath, therefore, is like a resonant 

cavity that Oscillates undetr..oiB):taln perturbations. .... 

4.2 Observatioiv-et-Plasma OscillalionS 

Initially, we looked for plasma oscillatiohs in the thin-sheath limit. An 
ambient density Of 2D0 cm*^ and an electron-ion temperature T^ * T^ = 0. 2 eV 
was chosen. This gives a Cebye length Of 0. 235 m compared with aim probe 
radius. The probe potential was biased at +1. 0 volt. The top three curves in 
Figure 10 show the sheath potential profile at 10 MseC intervals. Some time- 
dependence is observed but its coherent property iS not clear. The periodic 
nature of the time -dependence is enhanced by plotting the potential at a constant 
distance (r = 2. 60 m) from the probe surface as Shown in the bottom part Of 
Figure 10. The local potential Oscillates at a fi equency comparable with the 
plasma frequency (upg = 7. 99 x 10® rad/sec). The nature of this oscillation is 
further delineated by performing a computer experiment. 

4.3 \ Computer Experiment 

Figure 11 represents a computer experiment in which all the parameters 
except the ambient density remained constant. The ambient density in each curve 
is higher by a factor of 2 compared With the curve Immediately above it. The 
points represent 16 iteration averages and the error bai s the rms deviation from 
this average. The observation point was taken at approximately 16 Xp b Deby# 
length) in each run. The average and Standard deviation of the oscillation period, 
as estimated from these pUv ilSO given. In each instance the observed 
averaged period, T, is Shot‘< u n the plaSma period t * i. H x 10"^/ n sec). 

In order to estimate the oscillation Wavelength, X, the Ohe-dimensional plaSma 
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Figure il. Simulated Flasm a Ob clllattoiie. £ act! curve 
rdji^esehts the vdltage itfAe^hlstory for vaf*ldU 6 densitifee. 
NE uAd are the nUmbUr of coAi]$uier electrons and Ions 
reapdctlvdiy. 7 ithe Is lA micrdseconds. ^ denotes wave 
length which equals 2 */k 





dlBpo^aldrt relation” l8 need. Although this Is derived Tor a cnrteolnn geometry, 
we aoBume It also holds for the epherleal ease In the thin sheath limit. The rela^ 

tlon Is 


uip ^ + 3 < kT> k^/tn 


(3) 


>''here , ^ 

\ = 2it Ik. Q = Observed frequency, ° plas*"® frequency 

( kT) * mean electron energy 

m “ electron mass 

k = wavenumber . 

The reepecllve wavelehgttl., .« obtelned from Eq. (3). are aha«n on the rlgh.hand 
eiae of Figure ll. The oentral ..luea of T ».re uaeO. Beeults Imply . ...elengU. 

Oh the order of the eheath Ihlokned. (Bb ■ 0 ^). Since the ''"'rand’t 

at +1. 0 «oU and la close to lero at Bb- the analogy betneen the eheath and 

resonant cavity IS seen to be very cloSe. 

Plasma oscillations ate LandaU-dampfed by the transfer of wave energy o 

electrons traveling allghtly belo» the «av6 phase velocity, V • u / , “ ' 

present code »e artificially enhmiced the high energy tall of the J 

iron distribution bUl relalnsd the ooadltlon 3f/3V < 0. Under these condlllon. me 
»ave amplitude showu In Figure 11 significantly decreased indlcatmg. as expecle . 

enhanced damping. 

The osclUallons shovm in Figure U approach the long 
kXjj <1. in that case, there Is ah analytic expression 

It Is 


” for the rate of damping. 


(hV' 


EXP 


2k^ 


(4) 


Where 7 = decay rate (sec-^). u«p = electron plasma frequency. While there 
substantial flncluatloft th v ttom its sensitivity to k. the shortest decay time is 
190 gsec for n = BOO cm*^ run. This Is sufficiently long so that no significant 

damping Is expected to be seen in Figure 11. u ♦ * puan thp 

At long wavelengths the ^avfe phase velocity. Vp. Is much greater than the 

thermal electron velocity so that there are few electrons In the resonant reg . 
At shorter wavelengths (k^D ^ ” »“e wave phase velocity approaches the thermal 
ItTc velocity. T1» prenenc. of mnny eleclron. In tbl. velocity region repl^ 
dompe the plaemn «m,e. While the nbove reeuUn Imply the uoefulne.o of pie. 
etmltton technlquen In looking et colleotlve behavior, eomc centlon ehoold 

noted. 
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The basic pt‘oblem is that the observed oscillations cart be either physically 
stimulated or artifieally dHvOn by high frequency grid or computer •'rtoisC" that 
couples to the lower frequency plasma oscillation by aliasing. The computer 
code should be sufficiently free from noise so ^hat observed oscillations are con- 
sistent with analytic criteria. 3 Tlie most common technique for reducing noise is 
to treat the computer particles as having a finite Sizers and to periodically 
smooth the distribution function (n both velocity and spatial coordinates. 20 These 
Sophistications should delineate the origin and nature of the oselllatiohs. 


5. SI MMAIIY A^U-t:O^(:l.L’SI0^S 


Agreement has been shown between a time-dependent simulation code and 
steady-state solutions to the Vlasov- PoisSon equations. Transient voltage rise- 
times due to photoemlssiort were presented and Space-charge limiting effects dis- 
cussed. Finally, a computer "experiment" was presented that showed the presence 
of plasma waves. Caution is noted for possible aliasing effects. 

It IS concluded, based on the above results, that a simulation approach is a 
valuable and versatile method for dealing With complex. plaSma-related Space 
phenomena. In particular, additional features can be added aS building blocks 
with little modification of the existing code. Care should be taken, however, in 
distinguishing between Computer-related and physically-related effects. 
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A(^>endix A 


In this appeniUk the method used to simulate a Maxwellian photosheath is 
briefly discussed. The computation time is much shorter if the speed and direction 
of the emitted particles are determined rather than the Individual velocity compon- 
ents. The probability, p, that the speed, U, is less than or equal is given by 

p V -JT /° >^ ™ . „„ 

^ O 


where U and are normalised to the thermal velocity. J kT/m. . 

One of us (GKY) has inverted Eq. A 1 to give as. a function Of p over the 
domain 0< p< 1, .This domain was dlvided-lnto six. intervals and six empirical.. 
appt^ximations determined. For p - T. the approach of Hastings^ has been 
adopted. 

The semiconvergent series for p as p - l~ is 


p = 1 




or 


^o^ = -2 In (1 - P) + In (uj) + In | 4 in 



(A2) 


substituting the lead term In the second term and ignoring terms of order U. 

6 

^0 a -p) 4 In t -2 in a - p)] 4 in I (A 3 ) 


we have used 


t ^ -2 in (1 - p) 4 S| 4 aj l)i [ -2 In (l - p)] 


(A4) 


° 1 + a,yP 4 Bgp^ 4 BgP^ 


(A5) 


and 


For p -• 0 , the form can be derived formally 



(A6) 


Let 


i- 

X. 


> 1 .'' 

k- ■ 

f. 

I: 

r 


U = S + C- S® + C.S® . 

O 3 3 

Substituting Eq. (A7) iii Eq. (A 6) gives 


(A7) 


P = 





(Cj + C3 


2 





The Second and subsequent coefficients can be set to zero by appropriate choice 
Of the (that is, Cg = 0. 1 , Cg = 0. 0221428571, etc, ) 

Thus p ^ or S = (4.5ff)^/® p^^® . (A8) 


For p — 0^, substitute Eq, (A 8) ih Eq. (A7). In some of the intermediate Intervals 
of p, slightly modified analytic forms provided greater accuracy. 

Experiments show that photoelectronS produced inside the material surface 
are emitted with a cosine distribution relative to the surface normal. ^ A second 
random number, q, generates the appropriate angle. 

The numbers p and q are generated by the function, RANP, resident OH the 
AFGL CDC 6600. This generates random numbers uniformly between 0 and 1, 
excluding the end points. This function Is a multiplicative congruential generator. 
Our use of the random numbers is not affected by the limitations pointed Out by 
Marsaglia. ^ The "seed" for RAJlF-lS the quaSirandom bit string which is gene- 
rated by the computer's real time Clock. 
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